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Abstract 

It is often true that an invariant mass constructed from visible decay prod- 
ucts of a heavy particle may attain a maximum(or minimum) for a certain 
kinematic configuration only — this fact can be used to reconstruct relevant 
particle masses from observed decay product momenta of events near the in- 
variant mass endpoint. MSSM neutralino and chargino mass reconstruction at 
the LHC from multi-lepton endstates is illustrated by way of example. 
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1 Introduction 



Within the first few years of LHC data collection and analysis, a key issue will 
not only be the search for general signs of New Physics (NP) beyond the Stan- 
dard Model (SM), but also quantitative measurement of any NP particles produced. 
NP mass spectra, in particular, can offer an important handle on discriminating 
between different NP models. In the highly theoretically-motivated Minimal Super- 
symmetric(SUSY) Standard Model (MSSM), for example, many of the 100+ free 
input parameters, among which various relations are predicted by specific models of 
SUSY-breaking, are directly coupled to values of SUSY masses. In contrast to the 
situation at a lepton collider, where the center-of-mass (CM) collision energy can be 
precisely tuned to sweep through NP mass resonances, the LHC produces partonic 
collisions whose CM energies vary unpredictably from event to event, washing out 
potential resonance structures. Moreover, NP states may decay only partially to 
visible and detectable particles, as in the R-parity-conserving MSSM, where sparti- 
cles typically cascade down to SM states plus an even number of invisible Lightest 
SUSY Particles (LSPs). Any MSSM mass reconstruction technique applied to LHC 
data cannot therefore depend on the precise CM energy or total visibility of decay 
products. 

In answer to these demands, the phenomenological community has innovated a 
number of different mass reconstruction methods, the most standard and well-tested 
among these relying on measuring endpoints of various invariant mass distributions 
constructed from visible final leptonic and/or jet 4-momenta [lj. In the MSSM, for 
example, 3-body decays of the second neutralino to a pair of leptons plus the LSP, 
X2 — » P- + ^~x\) gives rise to a dilepton invariant mass distribution Mg+g- which cuts 
off relatively sharply at M'p^L = m^o — m^o . With a sufficiently large event sample 
this mass difference can be determined very precisely (to the sub-GeV level), but the 
values of the individual masses themselves are undetermined — this is typical of the 
endpoint method, where the endpoint is generally some function of a number of NP 
masses. 

The present work introduces the idea that there can be more information at the 
endpoint than just its numerical value — events lying at the endpoint may arise 
from a unique kinematic configuration of final momenta in a decaying particle's 
frame — and additional analysis can be done to find the masses. Let this be defined 
as the DK (Decay-frame Kinematics) technique. As a first example, consider the 
3-body decay noted abovcl An event at the dilepton endpoint must be such that, 
in the x° decay frame, the Xi is a t re st while the leptons are produced with equal 
and opposite momenta (the dilepton system has zero velocity). This fact allows 
us to find the velocity of the X2 m the detector frame: namely, it is the observed 
velocity of the dilepton system, (p 1+ + Pi-)/ (Eg+ + En- ) . When we apply the 
corresponding Lorentz boost to the LSP 4- vector (m^o, 0, 0, 0) and match to observed 
missing transverse momenta, the LSP mass can be easily solved for - - twice in 
fact, from each transverse direction. The fact that R-parity requires a pair of such 



2 The case of 2-body decays is quite similar and will be reserved for a future work [3J. 



decays in each event changes nothing essential: one simply finds both x® velocities 
from events at a double- endpoint, i.e. where both dilepton invariants are maximal, 
matching missing momentum to the sum of LSP momenta. Although no event ever 
lies exactly at a double-endpoint, a large number of events may be within a tolerably 
small neighborhood of it — one would expect that using these to reconstuct the LSP 
mass as above would give a distribution peaked at the correct value. There are, of 
course, immediate practical limitations — how small does the neighborhood have to 
be, do detector effects and backgrounds smear the peak beyond recognition, etc. - 
which must ultimately fall to a Monte Carlo test in order to properly address. 

As a second example of the DK technique in this paper, we shall consider pro- 
duction of neutralino-chargino pairs x^Xi m the MSSM, with subsequent decays 
X2 — * P + P~Xi an d Xi ~^ P ±l/ Xv Though slightly more complicated than the case 
of neutralino-pair decays, the same principle works: events at a particular endpoint 
belong to a certain class of kinematics in the X2 an d Xi decay frames, allowing us 
to find the individual velocities of these latter; the LSP 4-momenta can then be 
boosted and matched to observed missing momentum for only one (correct) value of 
the chargino mass. 

The structure of this work will be as follows: Section [2] will illustrate the DK 
technique in more detail in the case of neutralino-pair production, where the goal is 
LSP mass-reconstruction — this includes a Monte Carlo simulation to test how well 
the technique might actually work with LHC data. Section [3] then does the same for 
neutralino-chargino pair production, again seeing how well the chargino mass can be 
reconstructed. Section H] summarizes these results, comments on the relation to other 
mass reconstruction techniques in the literature, and discusses general applicability. 

2 Neutralino-Neutralino Modes 
2.1 Theory 

Consider production of neutralino pairs in the MSSM which undergo 3-body decays 
to electrons and muons: 

pp _> X _> X ' + e + e~x?) P Vx?) (1) 

where X represents either Z* or any MSSM production channel from a Higgs (H° or 
A ) or colored gluino/squark cascades, while X' are SM states potentially produced in 
association, all irrelevant to the current discussion; note we could include e + e~e + e~ 
and endstates as well, but at the price of a two-fold ambiguity in lepton 

pairing in what follows. Physical observables from one event thus consist of four 
leptonic 4- momenta p e ± M ± (from which we may construct the usual dilepton invariant 
masses, M ee and M w ) and missing momentum in two transverse directions, assumed 
equal to the sum of the two x? s ' transverse momenta, pf 2 - Decay kinematics (see 
Appendix) allow us to write the following list of constraints on the relevant neutralino 



masses (hereafter we abbreviate m, 

1 



Pe++ Pe- 



~P»+ + 



2m, 
1 

2m, 



mt 



m 2 ) 2 



4mfm 



(Mf 



mt 



m 2 ) 2 



4mfm 2 



~Pi + "p*2) T = (observed) 



(2) 

(3) 
(4) 



where leptonic momenta are written in the frame of the respective parent neutralino, 
i.e. p' e ± = Axp e ± and p' ± = A 2 p M ±, defining the appropriate Lorentz transformations 

A li2 . In general, of course, we do not know the Lorentz boosts (3 i j2 from which 
A li2 are constructed, so ©-(rjj represents a system of four equations (or six, if the 

dilepton endpoints M^ ax = m^j — m x are known) for nine unknowns ( (3 i )2 and the 
masses thi^j) which obviously cannot be solved uniquely for the masses. 




Figure 1: An event with maximal M ee and M^; though the decaying neutralinos • may 
be moving with any velocity {3\^ in the lab frame, in each respective decay frame the leptons 
have equal and opposite momenta while the x? is at rest. 



Supposing, however, we have an event where the invariant masses M ee and M : 
are maximal, as in Fig. [H an enlarged system of constraints results: 
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which now gives ten equations (([7]) and (IE]) stand for three constraints each) for the 
nine unknowns^], allowing us to actually overconstrain the masses The (3 i i2 

which satisfy ((?]) and (E]), making the total momentum of each lepton pair zero, are 
uniquely given by 
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3 Strictly speaking, there is a two-fold ambiguity between © and ©, unless i = j. 



Now the corresponding A 12 which take the e + e Xi an d A 4 4 I 1 Xi systems to their 
respective -rest frames also bring each xi to rest (a condition of M ee and M w 

being maximal): their 4-momenta in these frames must thus be (mi, ), which, when 

inverse-Lorentz-transformed by Aj~3, to give (m 1 7 12 , mi( (3 7)1,2)) must agree with 
the observed missing momentum j? T ; a matching condition along each transverse 
direction (say x and y) then gives two independent determinations of m\\ 

m > = m " - ~$y (\\\ 

Since we are assuming that both M ee and M MAt are precisely maximal (the perfect 
event of Fig. [T]), we would of course get m[ = m'[ = mi; m practice, of course, we 
can only expect to capture an event within some neighborhood e of the endpoints, 
M eetfll j, = M™™ ± e, in which case one can show (see Appendix) that m[ and m'{ 
are approximately ~ mi± \f2em\. One might then expect that applying ( JTUl) and 
fTTXH to a sample of events near the endpoint should give a distribution of m[ and m'[ 
peaked near mi with a 0(v / 2emi) spread. 

As a technical caveat to the above, things could go wrong if ~$ T = : this oc- 
curs if the neutralinos Xij happen to travel in opposite transverse directions with 
precisely the correct velocities, i.e. (^ x ,yl)i + {Px,yl)i — 0- In this case the formu- 
lae (Hip give us indeterminate solutions m[ = m" = 0/0. Even if is finite but 
small (say a few GeV) we might still worry that statistical fluctuations to ~ 
would likewise throw off the solutions pip . This will not be a serious concern, how- 
ever, if most events contain significant "upstream transverse momentum" (UTM) [3] 
in X or X', and we will presently see that a Monte Carlo simulation of a typical 
LHC environment, adding in experimental effects such as measurement errors and 
inherent finiteness of detector resolution, as well as the inclusion of SM and MSSM 
backgrounds, basically confirms the robustness of DK mass reconstruction. 



2.2 Monte Carlo Test 

In order to see how well the above programme might work with real data, let us 
apply it to Monte Carlo simulated LHC data, choosing for definiteness the following 
MSSM parameter point: 

H = 390 GeV tan(3 = 10 m A = 400 GeV 
Mi = 100 GeV M 2 = 123 GeV m- g = 605 GeV 
m~„ = 500 GeV m T . = 300 GeV m r = 130 GeV 

which is a modification of the SPSla' benchmark point [1] enhancing the x.2 ~* 
£ + £~Xi branching ratio. 

Colored sparticles, dominating the bulk of the inclusive SUSY cross section (~ 
70 pb) at this point, cascade predominantly to X2 (~ 30%) or xf (~ 60%), excepting 
R-handed squarks (which mostly decay directly to the LSP). As can be seen from 
the gaugino mass spectrum in Table [TJ the mass splitting between %° and Xi (the 



Table 1: Gaugino masses (in GeV) at the MSSM parameter point considered. 
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LSP) of m,2 — 777i = 18.7 GeV is too small to allow anything but xl decay through an 
off-shell intermediary (slepton or Z*) as desired: xl ~* ^ + ^~Xi nas BR ~ 90% (next 
of importance is xl ~^ Ql'Xi with BR « 9%). SUSY e + e~/i + /i~ endstates should 
therefore arise almost entirely from X — > X' + e + e~x?) XiK - * / i+ A i_ Xi)- 

Proton-proton collisions corresponding to 100 of LHC luminosity are then 
generated with HERWIG 6.5 [5] (coupled to the CTEQ6 parton distribution func- 
tions [6]) for all relevant SUSY processes (pp — > any pair of {q, ^,x ± )X°}) — * ^0 
and SM backgrounds (for a hard 4-lepton signal only Z^*'Z is sizeable [7]), cou- 
pled to a simplified detector simulatocl which only passes events with four hard, 
isolated^! leptons with flavor structure e + e~fi + fi~. The wedgebox plot of Fig. [2b, 
shows all surviving events (~ 5000 of these) by position in (M ee ,M AtM )-space and is 
well-suited to picking out those arising from the desired x°x°-origin [TU] ; in this case 
the dense box-like structure in the corner (M ee>flfl < 20 GeV) arises from x^^-pair 
decays as expected; areas of the wedgebox plot not relevant to the present study 
include the 'wings' projecting along either axis, arising from decays of more sparsely- 
produced X2X3 4 an d H pairs, as well as the two lines of points concentrated around 
M eew » 90 GeV due to the Z&Z background. 

Focusing on the corner region then, a flavor-subtracted dilepton invariant mass 
distribution (Fig. [2b) constructed from hard 2-lepton events (over 27000 of these) is 
most useful for pinpointing the endpoint: for our purposes it is enough to roughly 
identify M^X a ^L « 19 GeV by visual inspection (this is already in excellent agreement 
with the expected value of m 2 — m\ = 18. 7 GeV). Note the less-than-triangular 
shape of this distribution does indicate the 3-body nature of the ^"decays responsi- 
ble. The X2X2 box on the wedgebox plot is of course defined where M ee>flfl < M]+^ , 
and an enlarged view of this region (Fig. [2b) allows us to select an event sample 
in the neighborhood of M ee)Wl ~ M™ a e _ — here 'neighborhood' might be defined as 
M™ a ^L — e < M ee)Wl < M^L for some optimal e, i.e. e must be large enough to get 
more than a few events, yet not so large that the kinematic configuration of Fig. [1] 
becomes a bad approximation for most events; at the present point e = 2 GeV seems 
best (see below), giving an event sample of size O(10 2 ). For each such event, ffTUl) 
and (TTTj) are applied from the measured leptonic momenta and missing transverse 
momenta, both m' and m" being added to the final distribution (Fig. [2bl) only if 
they agree with each other to within 20% — this is very effective at eliminating 

4 An identical set-up was employed in the author's previous publications [51 [S] and the reader is 
referred there for additional details. 

5 No tracks of other charged particles are present in a r = 0.3 rad cone around the lepton, with 
less than 3 GeV of energy deposited into the electromagnetic calorimeter for 0.05 rad < r < 0.3 rad, 
and p l T > 10, 8 GeV for I = e, //. 
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Figure 2: (a) Wedgebox plot (100 fb^ 1 ) showing a very dense box with diffuse backgrounds 
(except near the Z-pole). (b) Flavor- subtracted dilepton events give a sharp distribution with 
an edge of approximately Ml+VL « 19 GeV. (c) Enlarged view of the wedgebox within this 
edge; events for analysis are taken from the small boxed region of area 2x2 GeV 2 shown, (d) 
Distribution of reconstructed LSP mass, with Gaussian fit overlaid: m^o = 89. 9± 19.1 GeV. 

not only background events which by chance might fall in the sampling region, but 
also legitimate X2X2 events for which the DK technique would fail (mostly because of 
momentum mis- measurement effects in the detectors but also possibly from wrongly- 
assumed kinematics of Fig. [1]) — nevertheless, some events give consistent solutions 
which lie far from the prominent peak at 80 — 100 GeV. At any rate, a Gaussian fit^to 
the peak indicates a reconstructed LSP mass of mi = 89.9 ± 19.1 GeV, in safe agree- 
ment with the nominal value in Table dJ The slight downwards bias in the central 
value seems to be a systematic effect, for increasing e enhances the bias — indeed, 
at lower e (selecting a sample of events more closely approximating the kinematics of 
Fig. [1]) the central value more closely approaches the nominal one, but at the price of 



6 For lack of a known theoretical shape. Away from the peak the distribution is clearly non- 
Gaussian, a feature which the author has found true at other parameter points tested. 



statistics, the balancing of these two effects leading to the quoted value of e = 2 GeV. 
Recall from the theoretical discussion above that the expected width of the LSP-mass 
distribution would then be at least 0{y/2^T l ) = 0(^(2) (2) (90)) GeV » 19 GeV, in 
good agreement with what is seen; a higher luminosity sample (even at full design 
luminosity of 300 fb^ 1 ) might permit a smaller optimal e, say e ~ 1 GeV, but the 
distribution width of ^/2em 1 would not be expected to decrease significantly — appli- 
cation of DK therefore reliably locates the LSP mass, but only to a limited precision 
(probably 10% at best). This may, nevertheless, be entirely sufficient for discrimi- 
nating among different SUSY models, and, as we shall see next, is also good enough 
to find the mass of the lighter chargino. 



3 Neutralino-Chargino Modes 

The foregoing has demonstrated how the DK technique reconstructs the LSP mass 
from 4-lepton endstates of neutralino pair (X2X2) decays; let us now see how DK 
also reconstructs the chargino mass from 3-lepton endstates of neutralino-chargino 
(X2X1 ) decays, again taking these to be through off-shell intermediates, i.e. 

pp _> X - X' + e+e-x?) XT"(- ^Xi) (12) 

or with e <-> /i; again, same-flavor endstates such as e + e~e ± or /i + /i~/i ± could also 
be included in this analysis at the price of a two-fold ambiguity in lepton-pairing. 




Figure 3: Events of interest for reconstructing the chargino mass: in the frame where 
the X2 decays at rest, the e are along some global coordinates ±y with maximal energy; 
likewise in the xf decay frame the muon is along ±x with maximal energy. The velocities 
Pi,2 are a priori unknown. 

With the X2X2 modes of the last section, it was natural to classify 4-lepton events 
by their dilepton invariant masses M ee , w in the form of a wedgebox plot, events of 
interest being concentrated near M ee = = M$ ax . Here with 3-lepton events, 
a natural way to classify events is by the dilepton invariant M ee and the invariant 
M121, defined as 

M m = {(Pe++Pe~ -P M ±) 4 + (ZV± +Pe+ ~ Pe'f + (Pe~ + P^± ~ Pe+f} /3 (13) 

since, as found in [S] , events with maximal M ee = M^ ax and minimal M; 2 ; arise from 
events as shown in Fig. [31 i.e. if the decaying xt an d Xi na d no motion then the e ± 



would have equal and opposite momenta (along say ±y) while the muon would be 
emitted perpendicular to y, say along x, also with maximal kinetic energy. 

Just as before, the required Lorentz transformation Ai to get the e + e~ pair in 
the frame of the decaying neutralino (p ' ± = A\p e ±) is determined by 
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(14) 



There is unfortunately no such simple formula for (3 2, but it is nevertheless uniquely 
defined by the following system of equations describing conservation of the total 
missing 4- momentum j^, 
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of which the two transverse components ~p* T are measurable, in addition to the three 
kinematic constraints 
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which events of the form shown in Fig. [3] must obey. Thus, ffl5l) and ffl6|) compose 

a system of five equations for the six unknowns { f3 2, mi, m.2, m±\. Supposing m\ 
is already determined (say from the last section), then we are less one unknown and 
the chargino mass m± can found from one (perfect) event, in principle. The usual 
practical caveats apply to this claim, so again we will have to check via Monte Carlo. 



3.1 Monte Carlo Test 

Working with the same MSSM parameter point and Monte Carlo setup as in the last 
section, signal and background processes (now including significant 3-lepton sources 
such as W^*/Z) for 100 jb~ x integrated luminosity are generated and now filtered for 
events with only three hard and isolated leptons — this, in addition to the fact that 
squark BR's to charginos are typically twice those to neutralinos at this point, gives 
a much larger signal: the "Mg+^- versus M121" plot of Fig. Hh represents well over 
4 ■ 10 4 events. Though this plot gives us some discriminating power between signal 
and background, e.g. the majority of SM background events cluster about the line 
M1+1- = mz (cut off by the plot), there is still a more-or-less uniform distribution 
of background events such as X2X1 rzy X?) (the t au decaying leptonically) which 
will have to be tolerated. 

From the previous theoretical discussion we know where to look on this plot, viz. 
the small rectangular box in Fig. containing events assumed to be of the type 
shown in Fig. [3j by trial and error it is found that the several hundred events here 
satisfying M^+i- = 18 ± 1 GeV and M121 < 20 GeV give an optimal sample in the 
sense discussed above in Section [2J For each such event, mi is chosen in the interval 




Figure 4: (a) Correlation between M^+£- and M121 (400 fb~ 1 ); events are taken in the 
boxed region shown where M^+^~ = 18 ± 1 GeV and Mi% < 20 GeV. (b) Distribution of 
reconstructed chargino mass from events in this boxed region, assuming the reconstructed 
LSP mass of mx = 89.9 ± 19.1 GeV. A Gaussian fit gives m~± = 106.2 ± 11.6 GeV. 



rrii — 89.9 ± 19.1 GeV (flat priors) and the system of equations ()14j) - (jl6p is solved for 
m±. Only physical solutions (m± is real and > mi) are kept and plotted in Fig. HJd. 
A Gaussian fit to the peak gives m± = 106.2 ± 11.6 GeV, in good agreement with 
the actual value m± = 114.2 GeV; note again, however, the downwards bias in the 
central value. As with the fit to the LSP mass in the last Section, repeating the 
analysis with various smaller-sized sampling regions decreased the bias at the cost of 
statistics, so this is understood to be a systematic uncertainty of DK. 



4 Discussion and Summary 

We have now seen two specific examples in the MSSM (from X2X2 an d xtxt P a h 
decays) of how DK operates: events near the endpoint of a particular invariant mass 
distribution (in the black boxes of Figs. Eb and Hk, respectively) can be used to 
reconstruct masses (mi and m±) to reasonable (~ 10 — 20%) precision, at least with 
the LHC luminosity and MSSM input parameters assumed here. Though 3-body 
decays xt ^ + ^~Xi an d Xi ~^ ^ ±l/ Xi were assumed for the sake of simplicity, 
the DK technique also works (with minor adjustments) for 2-body decays through 
on-shell sleptons [2]. 

It is important to ask how this technique performs relative to other mass re- 
construction methods available (excepting, of course, the traditional endpoint-value 
method, which is incorporated into DK). First consider the reconstruction of neu- 
tralino masses. There are now an array of methods which take advantage of the 
pair-production of neutralinos. One class of "Mass-Shell Techniques" (MST) , repre- 
sented in the work of [TT] and [T2] , essentially depends on maximizing the solvability 
of assumed mass-shell constraints in a given sample of events. This seems quite 



effective for on-shell decayCl but for the off-shell decay topologies in the present 
work these methods cannot be applied since there are not enough such constraints. 
Recently fashionable "transverse mass variable" methods [131 EL e -g- m T2, might 
be applied. It is, in fact, quite likely that such methods are closely related to DK in 
the sense that they also attach importance to events with particular kinematic con- 
figurations, e.g. the analysis in [2] would suggest that a trti variable constructed 
for — * e+e f i+ t I ^Xi attains a well-defined (in terms of sparticle masses) max- 
imum when the lepton pairs are either back-to-back or collinear and parallel; in one 
such development [3] , for example, a "constrained mass variable" m,2c in events with 
large UTM proves quite powerful for even a small number of events. In any case, if 
any of the above methods are applicable, they should perform much better than DK 
alone, since all events (not just those close to a double-endpoint) are used. As 
for decay modes with charginos, such as the XiX^, modes considered in this work, it 
would seem that other techniques would encounter difficulties due to extra invisible 
particles (neutrinos) in the decay products. There is the interesting prospect, how- 
ever, of somehow conjoining these techniques with DK for a better fit of the chargino 
mass. 




Figure 5: General 2d Dalitz plot; (a)-(d) indicate special regions which may have unique 
kinematics applicable to DK analysis. 

The DK technique is by no means limited to neutralino or chargino decays, nor 
just to MSSM applications. Any NP decay chain which is subject to the traditional 
invariant mass endpoint method also can be analyzed with DK. The key is to choose 
invariant masses whose endpoints correspond to unique decay-frame kinematics, e.g. 
M«(at its maximum) in the case of neutralino decays, or both Mu and M121 (which 
must be simultaneously maximal and minimal, respectively) for neutralino-chargino 
decays. This statement allows one to imagine how the DK technique might work 
in general: depending on the specific number and type of invariant masses one can 

7 But see [10] for some important caveats. 



construct from the endstates of a given decay chain, some d-dimensional Dalitz plot 
of these will be useful for isolating events with a unique decay-frame kinematic con- 
figuration. Fig.[5]shows, for example, a hypothetical 2d Dalitz plot of some invariants 
M a and M$: points which might exhibit unique, identifiable decay-frame kinematics, 
hence amenable to DK analysis, include (a)extrem£(§, (b)kinks, (c)internal density 
minima or (d) maxima. Extension to d > 3- dimensional Dalitz plots is feasible if 
enough events can be collected to fill in the shape. 
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But beware: endpoints do not always arise from unique kinematic configurations; see [51 115j. 



Appendix 



Off- Shell Kinematics 

Consider the three-body decay Xi ~ > ^ + ^~Xi m the rest frame of the decaying neu- 
tralino. Working in the approximation where the lepton mass is zero, the leptonic 
four-momenta are p^ ± = (E±, E±p±) and the dilepton invariant mass Mf +E _ = 
(E + + EL) 2 — (E + p + + E_p_) 2 can be rewritten as 

il^V = 2E+E_(1 - cos9) (17) 

where 6 is the angle between the leptons. Energy and momentum conservation further 
imply 

E + + E_ + Ei = rrii (18) 
E + p + + E_p_ + ~p l = (19) 

where the four vectors of the X12 i n this frame are (Ei, ~pi) and (m*, ), respec- 
tively Equations f|T7j) - f|T9l) give 



\E +P+ + e-p.\ = Vm^E^E^E^ (20 ) 

which in the limit where M^+i- attains its maximum (= — mi) reduces to 

\E+p+ + E-p-\ =0 (21) 

further implying 

rrij - mi , . 

E± = -^—^ (22) 

cosO = -1 (23) 

i.e. the one constraint (I20p is replaced by (j21j) - (l23j) for events with maximal dilepton 
invariant mass. 



Spread of LSP mass measurement 

In the discussion of Section [2] we showed that if we have the perfect event M ee = 
M w = MJl ax then the extracted masses in ffTTT) both give m' 1 = m'{ = mi. In a finite- 
size event sample we can only expect to approach perfection, say M ee)W , = M$ ax — e. 
In this case the Lorentz transforms defined in (fTUI) will not give the LSP rest frames, 
where by definition their energies are E\^ = mi, but rather frames where 



as can be derived from (I17p - (jl9p above, with corresponding momenta 
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Inverse-Lorentz-transforming these four-vectors gives 
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for each {3 = (3 \^ which modifies (II II) to 

~fx 



m. 



Pixli + 02x12 ± 0( v /2e/m 1 ) 



Ay 71 + P2yl2 ± O(^S) 



(27) 



The precise spread in the distribution of values of mi therefore depends on a con- 
volution of the Xi'-vslocity-distribution and three-body phase space, but from the 
above we can already see it is roughly of order y/2mi\. 



References 

[1] H. Bachacou, I. Hinchliffe, and F. E. Paige, "Measurements of masses in SUGRA 
models at LHC." Phys.Rev.D62:015009 (2000); E. Lytken, "Derivation of some 
kinematical formulas in SUSY decay chains." ATLAS note ATL-PHYS-COM- 
2004-001; J. M. Butterworth, J. Ellis and A. R. Raklev, "Reconstructing sparticle 
mass spectra using hadronic decays." JHEP 0705, 033 (2007); B.K. Gjelsten , 
D.J. Miller , and P. Osland. "Measurement of the gluino mass via cascade decays 
for SPS la." JHEP 0506:015 (2005); P. Huang, N. Kersting, and H.H. Yang, 
"Extracting MSSM Masses From Heavy Higgs Decays to Four Leptons at the 
LHC." Phys. Rev. D77, 075011 (2008). 

[2] N. Kersting, work in progress. 

[3] A. Barr, G. G. Ross, and M. Serna, "The Precision Determination of Invisible- 
Particle Masses at the LHC," Phys.Rev.D78, 056006 (2008). 

[4] B.C.Allanach et al., "The Snowmass points and slopes: Benchmarks for SUSY 
searches," Eur.Phys. J.C25, 113-123 (2002). 

[5] G. Corcella et al, JHEP 0101: 010 (2001); S. Moretti, K. Odagiri, P. Richardson, 
M.H. Seymour and B.R. Webber, JHEP 0204: 028 (2002). 



[6] J. Pumplin et al. (CTEQ Collaboration), JHEP 0207:012 (2002); 
D. Stump et al. (CTEQ Collaboration), JHEP 0310:046 (2003); 



[7] G. Bian et al, "Wedgebox analysis of four-lepton events from neutralino pair 
production at the LHC," Eur.Phys. J.C53, 429 (2008). 

[8] N. Kersting, "On Measuring Split-SUSY Gaugino Masses at the LHC," 
larXiv:0806.4238l [hep-ph] 

[9] P. Huang, N. Kersting and H.H. Yang, "Extracting MSSM masses from heavy 
Higgs boson decays to four leptons at the CERN LHC," Phys.Rev.D77, 075011 
(2008). 

[10] M. Bisset, N. Kersting, and R. Lu, "Improving SUSY Spectrum Determinations 
at the LHC with Wedgebox and Hidden Threshold Techniques," arX iv: 0806 .24921 
[hep-ph]. 

[11] K. Kawagoe, M. M. Nojiri and G. Polesello, "A new SUSY mass reconstruction 
method at the CERN LHC," Phys.Rev.D71, 035008 (2005). 

[12] H.-C. Cheng et al, "Mass determination in SUSY-like events with missing en- 
ergy," JHEP 0712, 076 (2007). 

[13] C. G. Lester and D. J. Summers, "Measuring masses of semi-invisibly decay- 
ing particles pair produced at hadron colliders," Phys.Lett.B463, 99 (1999); 
A. Barr, C. Lester and P. Stephens, "m(T2): The truth behind the glamour," J. 
Phys.G29, 2343 (2003). 

[14] W. S. Cho et al., "Measuring superparticle masses at hadron collider using the 
transverse mass kink," JHEP 0802, 035 (2008), Phys. Rev. Lett. 100, 171801 
(2008); A. J. Barr, B. Griparios and C. G. Lester, "Weighing WIMPS with kinks 
at colliders: invisible particle mass measurements from endpoints," JHEP 0802, 
014 (2008); M. M. Nojiri et al., "Inclusive transverse mass analysis for squark 
and gluino mass determination," JHEP 0806, 035 (2008). 

[15] B. K. Gjelsten, D. J. Miller, and P. Osland, "Measurement of SUSY Masses via 
Cascade Decays for SPS la," JHEP 0506, 015 (2005). 



